Phonon dispersion and low energy anomaly in CaCg. 
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We report measurements of plionon dispersion in CaCe using inelastic X-ray and neutron scatter- 
ing. We find good overall agreement, particularly in the 50 meV energy region, between experimental 
data and first-principles density-functional-theory calculations. However, on the longitudinal dis- 
persion along the (111) axis of the rhombohedral representation, we find an unexpected anti-crossing 
with an additional longitudinal mode, at about 11 meV. At a comparable energy, we observe also 
unexpected intensity on the in-plane direction. These results resolve the previous incorrect assign- 
ment of a longitudinal phonon mode to a transverse mode in the same energy range. By calculating 
the electron susceptibility from first principles we show that this longitudinal excitation is unlikely 
to be due to a plasmon and consequently can probably be due to defects or vacancies present in the 
sample. 
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I. INTRODUCTION 

Intercalation of foreign atoms in graphite can stabilize 
superconductivity by introducing metal atoms between 
the layers, which allow both tuning of the inter-layer 
spacing and charging of the graphite host. For a long 
time it was believed that the maximum critical temper- 
ature obtainable at ambient pressure in graphite inter- 
calation compound (GIC) [1] was less than 2 K. The 
discovery of high temperature superconductivity in two 
intercalated compounds: YbC6t2|] and CaCed, Q with 
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unprecedented high transition temperatures, 6.5 K and 
11.5 K respectively, has raised renewed interest about the 
role of phonons in GICs. 

The role of the in-plane and out-of-plane phonon 
modes in CaCg has been controversial. By using the den- 
sity functional theory, it has been shown that supercon- 
ductivity in CaCg is due to an electron-phonon mecha- 
nism y-ISl • The electron-phonon coupling is mainly asso- 
ciated to Carbon vibrations perpendicular to the planes 
of graphite (Cz), and Calcium vibrations parallel to the 
graphite planes (Caxj,)- However specific heat measure- 
ments suggest Q that the contribution of Cz vibrations 
is even larger than what was predicted by density func- 
tional theory (DFT) . Consequently, it would be desirable 
to measure the phonon dispersion of CaCe in order to see 
if there is actually a disagreement between theory and 
experimental data. 

In this work we measure phonons in CaCe along the 
(111) axis of the rhombohedral representation using in- 
elastic X-ray and neutron scattering (IXS and INS re- 
spectively). Furthermore, in order to asses the precision 
of DFT simulations we compare the calculated and mea- 
sured in-plane averaged IXS structure factor. 

On the longitudinal dispersion along the (111) axis of 
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FIG. 1: (Color online) Right panel: rhombohedral crystal 
structure of CaCe- Large red circles correspond to calcium 
ions, small black ones to carbon atoms. Left panel: The cor- 
responding hexagonal unit cell. The 3- rhombohedral c-axis 
is perpendicular to the honeycomb graphene planes. 

the rhombohedral representation, we find an unexpected 
anti-crossing with an unidentified additional mode, at 
about 11 meV. At the same energy, we observe also un- 
expected intensity on the in-plane direction. 

Previous measurements of phonon dispersion in 
CaCe [1| assign a longitudinal character to a mode in the 
same energy range, but the mode is transverse. Here we 
show that and additional mode appear in the same en- 
ergy region, and interact with the longitudinal acoustic 
mode. This result could help to resolve this apparently 
incorrect mode assignment. 



II. METHOD 

IXS and INS measurements have been conducted on 
a polycrystalline sample of highly oriented pyrolytic 
graphite intercalated with Ca, prepared by liquid-solid 
synthesis method in lithium-based alloys H • The sample 
shows same superconducting properties [lO| as previously 
reported on samples prepared using the same method [3| . 

Phonon measurements in graphite are a classical exam- 
ple of comparison between IXS and INS, as shown at the 
very beginning of the IXS technique by E. Burkel [Tl| . 
The IXS approach is particularly favorable for high en- 
ergy modes as it is the case in graphite and GICs where 
some optical branches involving carbon ions are in the 
170 meV energy range [12i] . 

In this kind of sample the polycrystalline domains are 
actually aligned along a crystal direction, corresponding 
to the the (111) axis in the rhombohedral structure of 
CaCe [ill (see FigH]). As a consequence, the phonon 
dispersion of the longitudinal modes, with momentum 
parallel to the (111) direction, can be measured. The 
different Carbon layers in the sample are randomly ro- 
tated respect to the common (111) axis so that the IXS 
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FIG. 2: (Color online) Top panel: section of the rhombohe- 
dral reciprocal space corresponding to the a*b* plane of the 
graphite sub- lattice in CaCe. Circle corresponds to fixed Q 
lines in the plane. Two F points in different hexagons are 
shown. The line connecting the two F points is parallel to 
(1 — 10). The large, dashed hexagon indicate the projection 
of the rhombohedral Brillouin Zone, with X the zone bound- 
ary along (1 -1 0). Bottom panel: the circles can be shifted 
in different Brillouin zones along the graphite c*, adding a 
component QzX(l 1 1). The L point lies halfway between 
two planes, at QzX(l 1 l) + (0.5 0.5 0.5). 



measurement for a given momentum of modulus Q in the 
a*b* plane will be an average over a circle of radius Q, 
as shown in Fig. [21 top panel. The volume of the sample 
we measured was enough for a measurement of the lon- 
gitudinal phonon dispersion along (111) using INS, but 
the signal from the in-plane average was way too weak 
for neutron, while strong enough for IXS. Therefore we 
decided to couple the two probes in order to achieve an 
extensive insight on the phonon structure of CaCg. 

The IXS experiment was carried out on the undulator 
ID28 beam-line at the ESRF. We have chosen to work 
with the Si (9 9 9) [H, [ll] reflection, with a wave-length 
of 0.6968 A-i (17794 eV) and an energy resolution AE 
= 3.0 ± 0.2 meV [l^. Additional spectra were collected 
using the Si (11 11 11) reflection, with a wave-length 
of 0.5701 (21747 eV) and an energy resolution AE 
— 1.5 ± 0.1 meV. The back-scattered beam is focused 
on the sample position by a gold-coated toroidal mirror, 
which provides a focal spot of h x v = 0.270 x 0.090 
mm FWHM. Further details of the used configuration 
are described in Ref. [l7|. The sample measured in IXS 
is a platelet with thickness of ^ 0.25 mm along (1 1 1) 
direction and an area of 4.2 x 0.9 mm. The sample was 
hold in a Lindemann glass capillary, sealed in a glove box. 

The inelastic neutron-scattering experiment was car- 
ried out on the cold source 4F2 and thermal source 2T 
triple-axis spectrometer at the Laboratoire Leon Bril- 
louin in Saclay, France. Harmonic contribution from the 
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FIG. 3: (Color online) Inelastic X-ray (top panel) and neutron 
(bottom panel) energy loss spectra atQ = G + q = 2x (1 , 1 
, 1) + (0.3 , 0.3 , 0.3), corresponding to a longitudinal phonon 
polarization along 0.6 of the V — L line. Data are fitted using 
a convolution of the instrumental function with an harmonic 
oscillator. In the top panel, theoretical calculations have been 
included (dashed lines). 



monochromator was reduced using a standard graphite 
filter, so final neutron wave- number was typically of 1.975 
A-i for 4F2, and of 2.662 A'^ or 4.100 A'^ for 2T al- 
though some scan at 1.975 was taken for comparison 
sake or for averting artifacts from monochromator har- 
monics. Collimation were of 60'-open-sample-open-open, 
and with a horizontally focusing monochromator and a 
vertically focusing analyzer. Several platelets, of about 6 
mm width and for a total of 2 mm thickness, were kept 
together in an aluminum foil. The sample was sealed in 
a aluminum can, sealed in a glove box, using an Indium 
gasket. 

A comparison of the IXS and INS data is shown for a 
point along the F — L line in Fig. [3] 

In parallel to IXS and INS data we carried out first- 
principles density functional theory calculations in the 
linear response flM [l9j . We use the generalized gradient 
approximation [20] and ultrasoft pseudo-potentials [2l| . 
Technical details are the same as in refs. 141. 1221 ■ 



III. EXPERIMENTS AND DISCUSSION 

Figure HI top panel, shows the measured phonon dis- 
persion along the c-axis using IXS and INS, as com- 
pared with first principles calculations. The low-energy 
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FIG. 4: (Color online) Top panel: Phonon dispersion along 
the (111) direction. Hollow and filled circles indicate, respec- 
tively, constant Q and constant energy INS scan. Triangles 
represent IXS constant Q scan. Lines represents ab-initio cal- 
culated phonon dispersion. The scans are of type QzX (1 1 
1) -I- (q q q) with and Qz = l,2,3, in rhombohedral lattice. Us- 
ing standard notation, we label LA and TA the longitudinal 
and transverse acoustic modes respectively. We label LO and 
TO longitudinal and transverse optical modes, respectively. 
Bottom panel: Zoom on the low energy part of the above dis- 
persion data. Circles indicate the experimental points (IXS 
and INS together). Dashed lines show the dispersion of the 
ab-initio longitudinal acoustic phonon and of the unknown 
mode. Continuous lines show the dispersion of the interacting 
mode described by Eq. ((2|, fitted to the experimental data. 
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FIG. 5: Real and Imaginary part of the susceptibility for 
CaCe. 



phonons dispersion along (1 1 1) is generally in agree- 
ment with theoretical calculations, with two notable ex- 
ceptions. First, the longitudinal optical mode is hardened 
in experiments (the hardening is 7 meV at zone border). 
Furthermore, although the sound velocity of the longi- 
tudinal acoustic mode is in good agreement with DFT 
calculations, above q = 0.1 and about 8 meV, the lon- 
gitudinal acoustic branch (LA) suddenly bends, and its 
energy at the zone boundary is 11.4 ± 0.1 meV, accord- 
ing to INS, which is very close to that of the calculated 
first transverse optical mode at the zone boundary L (la- 
beled TO in FiglH top panel), at 11.8 meV. The same 
behavior can be observed on the high energy section of 
the LA mode. Close to the zone boundary, from q = 0.5 
to q = 0.3, the dispersion, as measured with both INS 
and IXS, reproduces very well the calculated one. Then 
it departs from the theoretical prediction going towards 
the zone center, below q ~ 0.3 and 16 meV. Finally to- 
wards zone center the LA mode flattens to a value of 
« 10 meV, slightly lower than what observed at L. 

Note, that the calculated energy for the TO mode at 
zone center F is substantially higher, namely 13.9 meV. 

This behavior is reminiscent of an avoided crossing, 
or anti-crossing, where the longitudinal acoustic mode 
interacts with an unknown mode. To analyze the ex- 
perimental data, we suppose that the dispersion of the 
unknown mode is described by: 



6cos(27r(7) 



(1) 



where a and h are fitting parameters and (q q q), is the 
momentum in the rhombohedral coordinates. We con- 
sider a g-independent coupling A with the longitudinal 
acoustic phonon. Within this model, the excitation ener- 



gies, observed in the inelastic experiments, are the eigen- 
values of a 2x2 matrix: 



A hwi^Aiq) J 



(2) 



where hui^xiq) longitudinal acoustic phonon en- 

ergy, as calculated using DFT. By minimizing the mean 
square error between the computed and observed ener- 
gies, we obtain a good fit with a = 10.75 meV, b — —0.82 
meV and A — 2.06 meV, see Fig. |4] (bottom panel). 
As a consequence the detected anomaly can be explained 
by the coupling of the longitudinal acoustic phonon mode 
with a longitudinal excitation or defect mode of unknown 
origin. 

The first optical mode is transverse supposing R3m 
symmetry, and therefore silent in this configuration, as 
it can be seen in the simulated spectra in Fig. [3l An 
interpretation of the anticrossing could be that the afore- 
mentioned first optical mode would be activated by some 
symmetry breaking scattering, as from impurities, which 
would mix different polarizations. However, the disper- 
sion, flat within error bars, or even with a slight positive 
slope at zone center, is not compatible with the negative 
slope of the first optic. 

A second possibility, could be the coupling to a longi- 
tudinal low energy plasmon. Indeed in another layered 
material, MgB2, an acoustic plasmon mode was indeed 
detected at surprisingly low energy [1^. To address this 
issue we calculate the G = and G = bare suscepti- 
bility, namely 



Xq(^) = Xq(G = 0,G =0,C^) = 



1 



E 

k.nm 



|(V'kn|e^*'l''"|V'k4-qn')P(/kn - ./k+qm) 



ekn - Ek-fqm + UJ - 11] 



(3) 



where Q =:496.38 a?, is the unit cell volume with an = 
= 0.53 A is the Bohr radius, and /k„ is the Fermi func- 
tion for a band energy ekn- Since the calculation of Eq. 
|3]requires a very accurate sampling of the Brillouin zone, 
we calculate the matrix-element using first-principles cal- 
culations on a 6 X 6 X 6 grid and then we interpolate it all 
over the Brillouin zone using Wannier interpolation (23 - 
f26| . In order to obtain a good description of the CaCg 
first principles bands, we use 7 Wannier functions. Then 
we interpolate the matrix element and the band structure 
over a Nk = 150 x 150 x 300 k-points grid and calculate 
Xq(a;) at g = (0.25,0.25,0.25), close to where the anti- 
crossing occurs. The temperature in the Fermi function 
is T = 300K and the lorentzian smearing rj — 4meV. The 
results for the real (x ) and imaginary (x ) part of the 
susceptibility are shown in Fig. [5] 

A peak in the imaginary part of Xq(i^) would signal 
the occurrence of a plasmon excitation. In our case the 
imaginary part of the susceptibility is featureless. Thus 
no plasmon excitations are present in CaCg calculated 
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FIG. 6: (Color online) Constant-Q in-plane phonon density of states obtained using IXS (line with error bars) compared to 
ab-initio calculations (thick lines). The experimental elastic contribution is also shown (thin lines). The scans are obtained 
at G + q where G = Qz x (111) and Qz = l (top, left), z=2 (top, right), z=3 (bottom), while q lies in-plane and we give its 
modulus along the (1 -1 0) direction in reduced length units for each spectra. For both top panels the energy resolution is 3 
meV. For the bottom panels the energy resolution set-up is SmeV (left) and 1.5 meV (right). Calculations are normalized to 
the mode at about 50 meV in the left panel, while normalization is to the lowest energy mode in the right panel. 
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susceptibility. To cross-check our result we also consider 
the case of MgB2 where a plasmon was found [1^. In 
the MgB2 case we reproduce the occurrence of an acous- 
tic plasmon mode in agreement to Ref. '23'. As a conse- 
quence the 11 meV excitations is unlikely to be due to a 
longitudinal plasmon, but most likely to impurities. 

A possible origin of the impurity mode would be the 
presence of lithium in the preparation of the samples by 
liquid-solid synthesis method [9] . This hypothesis is sup- 
ported by a strong XPS signal for the lithium ion at the 
K-edge (results not shown). However, this measurements 
are not confirmed by nuclear microprobe analysis. 

A further possibility would be the presence of a 
PGa/mmc phase in the sample coexisting with the most 
stable R3m one. However the attempt to describe the 
c-axis phonon dispersion in terms of a similar P63 /mmc 
phase is also unsatisfactory since the dispersion of the 
low energy modes is hardly distinguishable from the dis- 
persion of the i?3m phase. Still it is possible that the 
avoided crossing is generated by the presence of vacan- 
cies or other similar defects. 

An anomaly at a similar energy is also observed in an- 
other set of data which consists in constant-Q 2D-phonon 
density-of-states, with a propagation vector q lying in 
the plane perpendicular to the (111) direction and cor- 
responding to the a*b* plane. For this IXS experiment, 
the choice of G = QzX (1,1,1), with ^ resulted 
in a significant contribution from transverse modes. In 
this respect what is measured is a 2-dimensional density 
of states obtained averaging over all the structure fac- 
tors having modulus of the IXS exchanged momentum 
|G -I- q|, where |q| is in the plane as reported in Figl2] 
and [6l Even if the measurement is not equivalent to a 
phonon dispersion calculation, as both longitudinal and 
transverse modes are measured, the low energy modes 
are fairly well reproduced, except for an extra intensity 
measured at about 10 to 13 meV for all Qz, in particu- 
lar for the spectra with q from 0.3 to 0.4, (see arrows in 
Figini in agreement to what we found along the c-axis. 

In a recent paper on a comparable sample of CaCg 
Q , using IXS only, Upton and co-workers incorrectly as- 
signed the flat band at '^ll meV to the first transverse 
optical mode. This is due to the fact that they lack both 
a symmetry analysis of their calculated modes as well 
as the resolution necessary to see the details of the anti- 
crossing. Note that the region close to zone center, below 
q = 0.15, is difficult to measure using IXS due to the large 
elastic peak at zero energy generated by disorder. This 
is not the case in INS scattering data and only by using 
this experimental technique the avoided crossing behav- 



ior between the LA mode and an unknown mode can be 
revealed. 

We note that for the first time, using the in-plane con- 
figuration with IXS, we found good agreement between 
the data and the simulation in the 50 meV energy re- 
gion. This is particularly important as vibrations propa- 
gating in-plane in this energy window are supposed to be 
among the ones contributing to electron-phonon mecha- 
nism P.[27i. 



IV. CONCLUSIONS 

In conclusion, we show that low energy phonon struc- 
ture of the superconducting graphite intercalated com- 
pound CaCe, show clear anomalies, in disagreement with 
DFT calculations assuming the perfect crystal structure. 
In particular we observe an anti-crossing of the acoustic 
longitudinal mode with an additional longitudinal flat 
mode. Calculation of the electron susceptibility suggests 
that this mode can not be attribute to a plasmon excita- 
tion. Thus we infer that the unknown excitation is most 
likely due to defects, impurities or vacancies. On the 
other hand, we found good overall agreement between 
the data and the simulations, particularly in the 50 meV 
energy region where the Cz vibrations propagating in- 
plane in this energy window are supposed to provide the 
largest contribution to electron-phonon mechanism. 
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